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I. Introduction 

Continuous phase modulation (CPM) is a widely used form of digital modulation which is em- 
ployed in many modern standards due to its very good power and spectral efficiency [1]. CPM is 
currently used in 2nd generation GSM (global system for mobile communications) cell phones, 
but also in IEEE 802.1 1 FHSS, Bluetooth, and other proprietary wireless modems. 
The (single-/^) CPM signal is defined by 

v{t) = e^'^W (1) 

where 

a(t) = 2-71 h'^anipit — riT) (la) 

n 

with: h the modulation index; {an} the data sequence with period T and M-ary alphabet Am = 
{±1, ±3, . . . , ±(M — 1)}; (^(t) the phase response, a monotonic function satisfying the property 
ip{t) = for t < and ip{t) = ^ for t > LT, with L the system memory. 

Thanks to its constant-envelope waveform, CPM guarantees a constant transmitted power, 
thus allowing use of low-cost nonlinear amplifiers, with a save of 3-10 dB loss of nonlinear 
distortion [1]. Moreover, due to the phase continuity, the CPM signals have good spectral prop- 
erties compared with memoryless modulations. The attractive feautures of CPM must face some 
major drawbacks, the most important being the high implementation complexity required for an 
optimal receiver, and the even more challenging task of receiver optimization in nonlinear and/or 
fading channels [2], such as satellite and mobile radio channels. A common approach to com- 
plexity reduction, proposed by Laurent in 1986 [3], is that of interpreting CPM as the cascade 
of an encoder followed by a bank of PAM (pulse amplitude modulation) interpolators, some of 
which carry very minor portions of the signal power and can thus be discarded. In this way, 
moderate-complexity efficient receiver can be built. Many different approaches of PAM de- 
composition have been proposed during years [4], [5], [6], [7], making CPM still an attractive 
research issue. 

More recently, the literature has been focusing on a generalization of the ordinary single-/i 
CPM format (1), formerly introduced by Anderson and Taylor in 1978 [8]. This generalization, 
known as multi-/i CPM, assumes that the modulation index is time-varying, to have 

a{t) = 27r^ an hni^it - nT) (2) 

n 
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where {hn} is usually a periodic sequence of period Nh. In comparison with single-/^ CPM, 
the multi-/i CPM signal is known to be able to reduce the bandwidth [10] and increase the 
probability of detection at the receiver side [1], [11]. Also, low-complexity receivers through 
PAM decomposition of multi-/i CPM signals have already been proposed in [7], [12], 

In this paper, we deal with another interesting aspect of multi-/i CPM signals, that is power 
spectral density (PSD) evaluation, an issue which has been investigated both in the past (e.g. see 
[9] and references therein) and in recent times [11], [13]. However, unlike the existing literature 
that has focused on numerical calculation, we aim at closed form PSD evaluation. 

This is possible by interpreting the CPM signal as a combination of PAM signals, using Lau- 
rent approach. The key point is then to model the encoder, preceeding the bank of PAM in- 
terpolators, as a finite-state sequential machine (SM) [14] whose output PSD can be derived in 
closed-form by exploiting the results of [15]. As will be clearer later on, the direct application 
of [15] is not possible in the multi-/i CPM context, since that theory is valid for a stationary 
and irreducible SM, whereas multi-/i CPM has a periodically time invariant (PTI) structure. 
Furthermore, the peculiar input alphabet Am leads to a reducible SM. In any case, a judicious 
modification of the SM case finally permits relying on [15]. Incidentally, the reducibility prob- 
lem is also found with the ordinary (single-/i) CPM case, of which we give new further insights. 
We also underline that the proposed method is closely related to the results of [16], with the non 
trivial differences of being able to propose a truly closed form and more compact results, and of 
being able to deal with a reducible SM. 

The paper is organized as follows. Section II deals with the representation of a multi-/i CPM 
modulator through a finite-state SM, while the investigation of PTI and reducible structures is 
given in Section III. Then, Section IV illustrates the formulation of a stationary and irreducible 
SM generating the multi-/i CPM signal after a proper bank of PAM interpolators. Section V 
recalls the findings of [15], and uses them to formulate the closed-form PSD evaluation. Some 
examples of spectra are reported in Section VI. Finally, Section VII concludes the paper. 
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II. Representation of CPM Modulator 

In this section we develop a convenient representation of a multi-/i CPM modulator, consisting 
of a SM followed by a bank of interpolators. The deduction of the CPM representation is based 
on the basic decomposition (BD) of the CPM signal into PAM waveforms [6]. The periodically 
varying modulation index leads to a PTI SM and to PTI interpolators. 

A. The basic decomposition for multi-h CPM 

The phase of the mvXii-h CPM signal (2) in the interval X„ = [nT, nT + T) can be written in 
the form 



a{t) = 271 



n—L 



amhrn ^ + <^mhm ^{t - niT) 

m=n— L+1 



m=—oo 



Then 



where 



L-l 



1=0 
n—L 



an-L= n With J^ = e^'^'^- 

m=— oo 

plays the role of a state, which is renewed according to the relation 



(3) 



(4) 



(5) 



jO-m + l 



(6) 



Comparison with the single-/i case [14] shows that the BD takes the form 



(7) 



where 



5. 



is a word sequence (column vectors of length = M^), and 

[L-l 



Quit) = vAt) 



Y[ exp {j27rhn-iai(f{t + iT)) 



i=0 



(7a) 



(7b) 



constitutes an interpolating filter bank (row vectors of length N). Here, ?7t(^) is an indi- 
cator function active over the interval [0,T), 5 is a vector generalization of the Kronecker 
delta function, a„ = (a„_L+i, . . . , a„) collects the input data over a window of length L, 
OL = (q!o, aL-i) G A\i, and is the cardinality of A\/i. 
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The difference with respect to the single-/^ CPM is in the renewal law (5), which is periodi- 
cally time-dependent (PTI), and in the pulses which are PTI with respect to the time n. 
This is shown in Fig. 1, where the non-linear device produces the symbols b„ from the input 
data and the bank of interpolators produces the PAM waveforms. 





PTI 




PTI 




an 


non 
linear 
device 




T Qui-) 


v{t) 


T 


T 


R 



Fig. 1 — A first model of a multi-/i CPM modulator. 

We now discuss the nature of the "state" am, whose renewal law is given by (6). We suppose 
that all the modulation indexes are rational and we write them in the form hi = Vi/p, where p 
is the least common denominator of the fractions. As an example [7], if /iq = \ = Jq^ hi = ^, 



h 



5 _ 10 



16' '"3 = 8" 16' '•^^ integer j9 is 16. Then we can see that 
It will be convenient to replace am by an equivalent state Zm such that 



which takes the values in the integer set 

e {0,l,...,2p-l} ^ U2p. 

The renewal law for Zm is obtained by rewriting (5) using (9), namely 



yj/2m gi'T^'m + iam + l/p _ +»'m + l Olm+l 



Hence 



^m+l — (^m + '^m+lO'm+l)2p 

where ( )2p denotes modulo 2p. 

In terms of Zm the word sequence becomes 



5n 



5a, 



n-L + l,aO 



,"L-1 



(9) 



(10) 



(12) 
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Now, the vector 6„ of length = M ^ can be conveniently written by the Kronecker product' 
®. Specifically, 

Proposition 1 Let Wa„ = [San ,a]aeAM be the indicator vector of a„ (a column vector of size 
M). Then 



For instance, for L = 1 and M = 4 we have 



(13) 



where 



tu_3 = [1 0]' 



w^i = [010 0]^ 



w 



+3 



[0 1]^ 



= [0 1 0] , 

B. Formulation of the sequential machine 

We now formalize the non-linear device of Fig. 1 as a PTI SM (see [14] for the single-/^ case), 
that is as a quintuple 

Mc,u = {A,B,S,g,h) (14) 

where A is the input alphabet, B is the output alphabet, S is the state alphabet, g is the state- 
update function and h is the output function, with 



(15) 



bn = h{sn,an,n] , B. 



Note that, in general, the state update function and the output function in (15) depend on n. 
When the dependence on n can be dropped we will talk of a stationary update/output function, 
and when the dependence is periodic on n we will talk of a PTI update/output function. 

^ Given two matrices A — 1 1 a,n„ 1 1 and B = 1 1 bij \ \ of arbitrary dimensions AIa x Na and Mb x Nb, the Kronecker product 
A X B is defined by 



A X B 



an B 



B 



and is a matrix of dimensions (MaMb) x (NaNb)- The Lth Kronecker power of a matrix A is defined asA^^ = AxAx 
■ ■ ■ X A (L factors). The following property relates the Kronecker product to the ordinary matrix product (mixed-product law) 
(A X B) (C X D) = (A C) X (B D) and more generally (Ai x ■ ■ ■ x A„) (Bi x ■ ■ ■ x B„) = (Ai Bi) x ■ ■ ■ x (A„ B„) . 
If A and B are invertible square matrices, (A x B)^^ — A^^ x B^^ . 
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We begin by expliciting the SM in a specific case, and let L = 3. The input is simply identified 
as a„ = a„, with alphabet A = Am. We then choose a valid state s„, and observe that am (or 
Zm) is not a sufficient information for determining (13). Hence, by inspection of (13), we let the 
state be 

«n = 4^^) = {Zn-3, an-2, On-l) (16) 

as in the single-/i case. Then, from (11), the state update function is expressed as 





= 




= a 


,(2) 


= a 



'2p 



,(2) 



(17) 



Therefore the state-transition has the form s„+i = g{sn, an, n) whose dependence on n is due 
to (17), where r„„2 depends on n. Moreover, the dependence is periodic in n because of the 
periodicity of r„. From (16) we see that the state set is 5 = N'2p x A?j^. Finally, the output 
function is defined by (13), which is explicitly 



w 



Wr 



JO) 



= W^2p "J^^Ci) ® W (2) ® Wa„ ■ 

This relation has the form bn = h (s„, a„), where h{-) is not time dependent. 

In the general case we have 
Proposition 2 In a multi-Zi CPM with alphabet size M, memory L and rational modulation 
indexes hn = Vn/p, n = 0, 1, . . . , A^^^^, the states of the SM Mcpu are defined by 



(1) o(i-l) 
1 ■ ■ ■ i^n 



(19) 



\Zn-Li O-n-L+li 

and the state set is 

s = ^^2px Air' ■ (20) 

The state-transition function Sn+i = g{sn, a-n, n), periodic in n with period Nh, is defined by 



(4°) + 



(1) 



2p 



(21) 



AL-^) _ JL-1) 
(L-1) _ 
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The output function bn = h a„) is given by 

^'n, = W^; ^W m ® ■ ■ ■ ® 1C (i^i) ^Wa„ . □ (22) 

The complexity of the SM TWcpm is essentially determined by the cardinality of the state set, 
given by 

I = \S\=2pM^-K (23) 
The input alphabet is .4 = v4m- The length of the output words b„ is 

l^^^l =M^ = N (24) 

and their structure is 

(0) 

[1 ••• 0] 

(0) 

[0 1 ••• 0] 

(0) 

[0 ••• 1] 
so that the output alphabet is a subset of >V^. 

C. Statistical description for spectral analysis 

The target is the evaluation of the average PSD Rv{f) of the multi-/i CPM signal v(t). In 
this analysis the following random processes are involved: the input data a„ = a{nT), the state 
sequence s„ = s{nT), the word sequence b„ = b(nT), and the CPM signal v{t), t G M. 

The assumptions on which we base our analysis are four, namely 

1) The input alphabet is an M-ary alphabet Am, with M even, containing odd symbols. 

2) The input data {a„} are stationary and statistically independent, with given probabilities = 

P [an = a],a e Am. 

3) The modulation indexes /i„ are rational, namely hn = rn/p,n = 0,1, Nh — l,p being the 
common denominator of the fractions and r„ being integers. 

4) The modulation index sequence {hn} is PTI of period Nh. 

As a consequence of 4), the sequential machine TWcpm of Proposition 2 and the filter bank of 
(7a) are PTI. This makes all the processes s„, and v{t) (Fig. 2). 
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Fig. 2 — Representation of a multi-/i CPM modulator by a 
PTI SM followed by a PTI interpolator filterbank. 



As we shall see, the period of cyclostationary depends not only on the period of the mod- 
ulation indexes Nh, but also on their parity. In fact, the cyclostationarity period is Tc = NcT 
with 

i Nh if Ef-n ^niseven 
Nc={ (25) 

[2Nh ifEn^O^niSOdd. 

Specifically, for a single-/;. CPM with even tq we have N^ = 1, but for a odd tq it is N^ = 2. 

The fundamental result for the statistical description of the above random processes is that, as 
a consequence of 4), the state sequence s„ is a non-homogeneus (i.e., non-stationary) Markov 
chain. Then, its full statistical specification can be obtained as a straightforward generalization 
of the results available for homogeneous Markov chains, providing the state absolute probabili- 
ties 

Pn{i) = P [s„ = i] , ieS (26a) 
and the state transition probabilities 

T^niiJ) = P [Sn+l = i\sn = j] , i,j eS . (26b) 

Both Pnii) and 7r„(i, j) are periodic of period Nc in n. They are customarily collected in matrix 
form under the name of, respectively, the state absolute probability vector (APV) p„ and the 
transition probability matrix (TPM) 7r„ defined as 



7Tr 



D. Evaluation of the transition probability matrix 

For the evaluation of the TPM we can use the same technique of [15] with the modification 
concerning the time-dependence. 
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6. 



(2^ 



Proposition 3 Let e^..^ = [ea,n{'>', j)]i,j£S be / x / binary matrices defined by the state- 
transition function s„+i = ^„(s„, a) according to 

ri if i = g^{j, a) 
lo if i^g^{i, a) 

We will call ea,n conditional transition matrices, since they express state transitions under the 
condition a„ = a. Then 

= J2 9» • ° (29) 

Now, by use of (28) and the state-update function (21), the conditional transition matrices can 



be explicited as 



^io,{jo+rn-L+lji)2p ^il,j2 ■ ■ ■ ^\ 



(30) 



io,ioeA/2p, il Jl,...,iL-l jL-lS^A 



where i = [io ■ ■ ■ and similarly for j. The product of functions in (30) permits formu- 

lating the matrix „ in compact form through a Kronecker product by following a procedure 
similar to that leading to (13). 

By preliminarily defining the single step ciclical shift matrix 

1" 



D 



2p 



1 



matrix 2p x 2p 



(31) 



which is a square matrix obtained by cyclically shifting the main diagonal to the left by 1 posi- 
tion, we then have (see Appendix A for a proof) 

Proposition 4 The I x I conditional transition matrices e^n for L > 2 are given by 



-■a,n ^2p 



As a consequence, and by exploiting (29), the TPMs 7r„ give 



w 



The case L = 1 must be treated separately. We have 



PeAM 



□ 



(32) 



(33) 



(34) 
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E. Evaluation of the absolute probability vector 

With respect to the APV the following update relation holds from the definition of p„ and 

Pn+l='^nPn- (35) 

The interest is to identify that particular APV which is invariant to the update (35), that is 

Pn = TTnPn , 

providing information on the limiting probabilities p^. Note that the APV we are looking for is 
an eigenvector of 7r„ with unit eigenvalue. As proved in Appendix B, such a vector of interest is 

1 , , 

P„ = — ® g ® ■ ■ ■ ® g • 36 

2p V ' 

L-l 



III. Classification of the Markov Chain 

A. Periodicity and reducibility properties 

We have already seen that the TPMs 7r„ are periodic with period A^^ . This makes the Markov 

chain being itself PTI. In addition, the Markov chain may further be reducible, that is if it is 

not possible to get to any state from any state or, equivalently, if some of the states cannot 

communicate between them. In the specific case of CPM, reducibility is assured by the rationale, 

explained in the following, and based upon the renewal law (11). 

Recalling that a,„+i is an odd number, we have the following prospect 

z„i even Zm+i even 

Zjn odd Zjn+1 odd 

z„i even Zm+i odd 

z,n odd Zm+1 even 

Now, for convenience in the state set S = x A'^J^ we partition A/2p in its even and odd 
parts 

U2r,{+) = {0, 2, . . . , 2p - 1} , U2r>{-) = {1, 3, . . . , 2p - 1} , (37) 

As a consequence, we arrange the matrix and vector indexes by displaying the even integers 
A/2p(+) first, followed by the odd integers A/'2p(— ), and find that the newly arranged matrices 
ea,n and 7r„ take a block diagonal form or a block antidiagonal form. 
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As an example, with L = l,p = 4, M = 2 and Vn-L+i = 3, we find tliat 7r„ has the block 
antidiagonal form 








2 4 


6 


1 


3 


5 


7 














9-1 


9i 





2 














9-1 


91 


4 








91 








9-1 


6 








9-1 


91 








1 





9-1 


91 










3 







9-1 










5 


q-i 


91 













7 





9-1 91 













e obtain the block diagonal fomi 











2 4 


6 


1 


3 


5 


7 








9-1 


91 










2 


91 


g-i 













4 





91 


9-1 










6 


9-1 


gi 













1 











9-1 





91 


3 








91 





9-1 





4 











91 





9-1 


7 








9-1 





91 






Note that submatrices are equal or closely related one each other by a simple circular shifts of 
the row (or column) order. 

The result can be generalized by inspection of (21) and (26b). Specifically, for tt^ we obtain 
the structure 

Tn-L+i even 

(38) 

Tn-L+l odd. 



7r„ = < 



F„ 



D'^J^Gr, 



where -F„ and G„ are /q x /q matrices with 



I, = -I = pM 



L-l 



(39) 



Note that all submatrices are themselves TPMs. 
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Moreover, submatices can be easily formulated as Kronecker products by recalling (33). For 
fn-L+i even we have 

with £)2p^^^'^(+l+) the square matrix collecting the samples of D2p'^^^'^^ at even rows and 
even columns, while for Tn-L+i odd it is 

with Z>2p~^'''^^(— 1+) the square matrix collecting the samples of Z)2p~^^^^ at odd rows and even 
columns. Note also that the correction term in (38) can be written as 

dJ/^ = Dp® I ML-. , (40) 

which restores the true modulo 2p operation when transiting from odd to even states. Inciden- 
tally, the following equivalences hold 

£>J/^ F„ = F„ D\f , D'jf G„ = Gn D\f . (41) 

By now decomposing the state set S = A/2p x into the subsets 

S{±) = U^A^) X Aij' . 
we can further illustrate the meaning of matrices in (38) by the following graph: 



rn-L+i even 



F 

s„ e — > s„e5(-) 



s„ e S{+) " s„ e S{-) 
r„_L+i odd D'^'^'G 

s„e5(-) : s„g5(+) 



(42) 



An equivalent illustration is given in Fig. 3. 

We see that, with Vn-L+i even, the state sets S{±) do not communicate: starting with s„ G 
S{+) the next state s„+i is again in 5(+), and this transition is governed by the TPM Fn- Anal- 
ogously, starting with s„ G S{—) the next state is s„+i G 5(— ) with TPM Fn- No transition is 
possible from 5(+) into S{—) and from S{—) into 5(+). For r.„_i+i odd, the state sets S{±.) 
do communicate deterministically: starting with s„ G S{+) the next state is s„_|_i G 5(— ) with 
TPM Gn, etc. 
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r„_L+i even r„_L+i odd 

Fig. 3 — State transition graphs from classes S{+) and S{—). 



B. Concluding remarks on Markov chain evaluation 

At this point we realize the fundamental role of the sequence r„_L+i of the normalized (to 
2p) modulation indexes. We start introducing two examples and conclude with the general case 
after. 

Example 1 Let L = 1, Nh = ^, tq even, ri and r2 odd. We consider the state sequence 
Sn Starting from n = 0. If sq G the transition is si G S{+) with TPM Fq, next the 

transitions are sa G S{-) with TPM d and S3 G S{+) with TPM D^J^' G2, etc. So we find 
the TPM sequence 

as illustrated at the top of Fig. 4. 



Q ro even j n odd 2 '"2 odd 3 ra even ^ r4 odd 5 rs odd g re even 




Fig. 4 — Trajectories of Example 1 . 

Analogously, if sq G 5(— ) we find the complemetary sequence illustrated below in Fig. 4. In 
both cases the period is N^. = = 3. 
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Example 2 Now let L = 1 and Nh = 3, and suppose tq and r2 even, ri odd. Depending 
on the state sq, we find the two trajectories illustrated in Fig. 5. In both case the period is 

= 2Nh = 6. 



ro even ^ n odd 2 3 '"3 even ^ n odd 5 rs even g re even 




Fig. 5 — Trajectories of Example2. 

Example 3 Let L = 1, Nc = 2, ro be even, and ri be odd. If sq G the sequence of 
TPMs is illustrated in Fig. 6 with period = 2Nh = 4. 



P ro even ^ n odd 2 '^2 even 3 rs odd 4 r4 even ^ rs odd g 




Fig. 6 — Trajectories of Examples. 

From the above examples and from the general rules (42) (see also Fig. 3), we see that, starting 
from n = 0, the Markov chain 

So , Si , S2 , . . . 

has two distinct classes of trajectories T(±) which do not communicate each other anytime. 
The first class T(+) is determined by the condition sq E S{+) and the second class T(— ) by 
the condition Sq G T(— ). We can extend these considerations to the bilateral chain 

which may be candidate for stationarity or cyclostationarity. 
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In the underlying probability space, we have to operate under one of the conditions^ 



C. 



and correspondingly we find two distinct classes of trajectories: 



{s„}gT(+) and {s„}eT(-). 



We claim that both classes can be modeled by a non-homogeneous irreducible Markov chain 
with TPMs 7r„(+) and 7r(— ), respectively. The cardinality of these conditional chains is pM^~^, 
that is half the cardinality 2pM^~^ of the original unconditional Markov chain. ^ 

The TPMs 7r„(±) are obtained from the TPM 7r„ of the unconditioned Markov chain with 
the block partition and the rules seen above and now summarized. 

Propositions Let 7r„ be the TPM of the unconditioned Markov chain, which has the form (38) 
and period Nc- The period N^. of 7r„(±) depends on the sequence of the normalized modulation 
indexes ro, ri, . . . , rAr,^_i, namely 



The 7r„(+) in a period starts with 7ro(+) = Fq if r„i+i is even (si G and with 7ro(+) = 

Go if r_2,+i is odd (si G S{—)). The rest is obtained recursively with the rule (42). Analogously 
7r„(— ) starts with 7ro(— ) = -Fo if r^L+i is even (si G and with 7ro(— ) = D^j°J^ Gq if 

r_i+i is odd (si G 5(+)), and is obtained with the same rules. □ 

IV. Alternative Representation of CPM Modulator 

A promising approach to closed form PSD evaluation is given by reinterpreting the CPM mod- 
ulator through the polyphase decomposition (PD) or serial-to-parallel conversion (S/P) of Fig. 7, 
where the input symbol sequence {a„} is turned into the word sequence 



which is itself a stationary sequence with alphabet A^f. 

■^In practice, the condition is determined by the remote history of the modulator. 

^The novelty (to be clarified) is that the state of each conditioned Markov chain is not unique but it changes from S{+) and 
5( — ) in dependence of the modulation indexes. However, T{+) and T( — ) are isomorphic. 




Nh if the number of the r„ odd in a period is even 
2Nh if the number of the r„ odd in a period Nh is odd. 



(43) 
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an 
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1 


Jn 


T <!>{■) 


v(t) 


T 


Tc 


Tc 


R 



PTI 



STI 



STI 



Fig. 7 

Representation of a multi-/i CPM modu- 
lator by a S/P conversion followed by a 
time invariant SM and interpolator filter- 
bank. 



The subsequent SM generating the output words can then be formulated on an equivalent 
to (3), with the only difference that now the observation interval of interest becomes of length 
Tc with the form X„ = [nTc, nTc + Tc). We then have 



a{t) = 271 



nNc-L 



nNc+Nc~l 



amhm 2 + o.mhm ^{t - mT) 

m=nNc — L+1 



m=—oo 



teXn 



(44) 



so that the counterpart to (7) becomes 

v{t) = Ycj,{t-nT)y^ 

n 

where 



(45) 



O'nNc-L 



L+Nc-l 



L-l 



Nc-1 



Yl exp {j27rhnNc-iai(p{t + iT)) J| exp {j27rhnNc+i/3i^{t - iT)) 



1=0 



i=l 



(45a) 



Here, rjTdt) is an indicator function active over the interval [0, Tc), an = {anN^-L+i, • • • , CLnN^+Nc-i] 
collects the input dataoverawindow of length A^c+-^—l5 CK = {ao, ai, . . . , ai-i, Pi, ■ ■ ■ /3nc-i) ^ 



Aij, and iVo = Af^+^'^~^ is the cardinality of ^ 



M 



So, is a row vector of length Nq, 



while 0(t) is a column vector of the same length. 

Note that, unlike (7), </)(t) is independent on n. Moreover, the second of (45a) clearly shows 
the newly required component driven by the symbols Pi. The formalization of a proper SM for 
(45) is now immediate by exploiting the results of Proposition2. 

Proposition 6 In a multi-/i CPM with alphabet size M, memory L, rational modulation 
indexes /i„ = r„/p, n = 0, 1, ... , Nh, and time-invariance period N^,, the states of the time 
invariant SM A^cpm (where p stands for parallel) are defined by 



r.= (af,ai^),...,af-)) = (z 



OnAfc-l) 



(46) 
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O-n = SnN, ■ (47) 

The output function y^ = h {cTn, Xn) is given by 

(0) 

Vn = ^2p W^^d) ® ■ ■ ■ ® W (L-i) ® W (0) ■■■(g) W ,N,-l) . □ (48) 

According to the ordering in (48) we can further attempt an equivalent Kronecker formulation 
of the interpolating filter bank cj){t) as 

0(i) = 0-L+i(i)®---®07v.-i(i) (49) 

where 

(t>i{t) = VtM J2 W^^hiatpit - iT)) (49a) 

In addition, by exploiting the sampling relationship (47), the statistical description of the 
newly introduced SM immediately follows from Propositions and Proposition 4, and from gen- 
eral properties of Markov chains. We have 

Proposition 7 The / x / conditional transition matrices Ec^^n for ct = [uq, . . . , aMc~i\^ and 
L >2 are defined through the matrix product 

(50) 

with „ as defined in (32). In turn, the TPM 11 is built as the TPM matrix product 

n = ''^nNc+Nc-l ■ ■ ■ '^nNc+l '^nNc (51) 

where 7r„ is given by (33). They both are independent on n. Observe that, for L = 1 expressions 
(34) must be used in place of (32) and (33). □ 

As a consequence of (38) and Propositions, we also have 
Proposition 8 By using the state ordering (37), the TPM matrix 11 takes the block diagonal 
form 



n 



n(+) 
n( 



(52) 
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where n(+) contains the samples of 11 at even rows and even columns, while n(— ) contains 
the samples of 11 at odd rows and odd columns. Moreover 

n(±) = 7r^,_i(±) ■ ■ ■ 7ri(±) 7ro(±) , n(+) = n(-) (53) 

the second equivalence being assured by the permutation property (41). □ 

Incidentally, an equivalent to Propositions holds for state-transition matrices E^, for which 
we have 

^-E„(+) 
E^{-)_ 

As a remark, we observe that the product formulation of (50), (51), and (53) could be restated 
through more direct expressions. However, this gives minor insights that are of no direct interest 
to spectral evaluation. 

Finally, the APVs of interest, separately for the even and odd trajectories T(+) and T(— ), 
can be formulated starting from (36) to give 

p = n(±)p, p = -lp®q®---®q (55) 

P ^ V ' 

L-1 

Incidentally, since both trajectories T(±) identify irreducible Markov chains, the following limit 
property is known to hold (e.g. see [17]) 

n(±) ^ \\mll{±f=pll, (56) 

that is the limit probability of being in a state is independent on the initial state. 



V. Closed Form Spectral Evaluation 

A. Generalities 

We follow the PD approach depicted in Fig. 7, which is more simple and permits the direct 
use of the theory of [15]. In fact, the reformulation of the SM in a time invariant form removes 
the PTI on the output sequence {y„} which, unlike {bn}, is stationary. Once evaluated the 
correlation 

ry{m = Y.\y^y*A (57) 
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and the corresponding PSD 

Ry{f)=T,T.ry{kn)e-''''^'^^ (58) 

k 

the output (average) PSD is simply given by 

RM) = ^if)Ryif)^*U-) (59) 

with $(/) = 4>{t)e~^'^'^^Ht the Fourier transform of cf){t). The result (59) is a straight- 
forward generalization of a widely known property of the interpolating filter to an interpolating 
filter bank, whose proof can be found in [ ???]. 

A fundamental role is played by the possible presence of lines (delta functions) in the SDs. In 
general, Ry{f) is given by the sum of a series, which is not convergent. The technique to handle 
this divergency is the separation of the correlation ry{kTc) into a continous part, r^y\kTc), and 
a discrete part, r'^\kTc). The latter is related to the asymptotic behavior of the correlation, 
and ultimately to the state transition probabilities (55). The PSDs will be henceforth explicitely 
written as 

Ry{f) = Ri'\f) + K'\f) (60) 

and 

Rv{f) = R^:\f) + Ri'\f). (61) 

We shall see that, when present, the spectral lines occurs at the frequencies multiple of Fc = 
B. New notation 

We preliminarily need to assess some notation according to [15]. Let the output symbols t/„ 
be organized in matrices such that collects (in its columns) the output words corresponding 
to the input word x. Evidently, being / the number of states, such words are in number of /, and 
matrix Yj, is thus a N x I matrix. By exploiting the first of (45a) and (48), in Kronecker form it 
is (see Appendix C) 

Ya, = V2p ® Iml-1 ® lu^o (g) ■■■ (g) w^^^_^ (62) 

where 

V2p = [HlieA/-^. = [1' ^2p, Wl, ■ ■ ■ , PFj"^] . (62a) 
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is a row vector. To discriminate between trajectories T(±), we can further introduce the notation 

Y^i+) = Vp (g) Iml-1 ® w^,, (g) ■■■ (g) w^^^_^ , Y^{-) = W2pY^{+) . (63) 

C. Presence of spectral lines in Ry 

Spectral lines depend on the limit behavior of the Markow chain, specifically they are related 
to the limit mean value of namely 

"^y = i™oE[t/J , (64) 
where APVs (55) hold. Now, the mean (64) can be evaluated as 

my{±) =J2F[x^ = x] Y^{±) p (65) 

X 

from which we have the following general result (the proof is reported in Appendix D) 
Proposition 9 When p > 1, no spectral lines occurr in CPM and we have 

i?W(/) = o, p>l. (66a) 

Conversely, when p = 1, that is for integer modulation factors, spectral lines are found and it is 

Ry^'\f) = ry^"^ E ^(/ - ^^c) , P=l (66b) 

A: 

where 



[qq*]®---0 [qq*] (67) 



andF, = l/T,. □ 

We note that Proposition 9 is in accordance to the results of the literature on single-/i CPM 
[18]. 

D. Continuous spectrum evaluation of Ry 

In order to evaluate the continuous part of the spectrum we preliminarily need to further 
explicit the correlation Vy in (57) following the development in [15]. 
For the value at nTc = (autocorrelation) we have (see (36) in [15]) 

ry(0) = E [y^ y*J=Y.^ [xm = x] YU±) diag(p) i;*(±) . (68) 

X 
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where, by exploiting the Kronecker product expressions (55) and (63), and after some little 
algebra, we have 

ry{0) = diag(g) ® ■ ■ ■ ® diag(g) . (69) 

V , y 

Nc+L-1 

For the values at n 7^ 0, we better distinguish between positive and negative values. For > 
we have (see (37)-(39) in [15]) 

ry{n%) = C2(±) U{±Y-^ Ci(±) , n > (70) 

where 



(70a) 



Ci(±) = 5] P K = x] E^{±) diag(p) Y^*{±) 

X 

C2(±) = EpK = ^]^-(±) . 

X 

Instead, forri < we can exploit the relation 

ry{-nTs) = ry\nTs) . (71) 

Observe that, the correlation samples of the two trajectories coincide since Ex{+) = Ex{—) 
and n(+) = n(-) (see (53) and (54)), while the W2p difference between Yx{+) and Yx{-) 
(see (63)) is removed by the complex conjugation in Ci. 

In addition, we must take into account for the limit value ry^'^\ which (when not null) is 
responsible of spectral lines. According to (40) in [15] we can write 

r^('^) = myi±)my*i±) = C,i±)UC,{±) . (72) 

As the reader can verify by exploiting (56), this expression is not in contrast with the formulation 
in (67). 

Now, the continuous spectrum becomes 

Ry^'^Hf) =TcY. (^y(^^c) - ry^'^) e'^^'^^"^^ (73) 

n 

and, by then relying on the symmetry in (71) and on the results of [15], the PSD can be written 

as 

i^/H/) = K(0) - ry^'^) 



(74) 

c,{±) {ij, - n) {e^'-f^^ij, - n(±) + n)~^ c,{±) 
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Note that this is a closed-form, unlike the results available in CPM literature (e.g., see [9], [13]) 
where the PSD is evaluated through direct numerical evaluation of the series (73). 

We also underline that the matrix inversion in (74) can be circumvented. The details can be 
found in [15], while here we report the final result. In particular, the inversion of an / x / matrix 
(AI — F), with A = e^'^^'f^", can be performed by 

(AI - F)-i = ^'=' ' , , G,= Y. dk-m 

"l^J m=0 

where 

d{x) = trace (xl — =^dk x^~^ . 

k=0 

E. Power spectral density of v {t) 

The derivation of the power spectral density of v{t) can now be straightforwardly obtained by 
use of (59). For ease of computational evaluation, here we introduce the matices 

K(/) = *(/)^.(+) (75) 

which are row matrices of dimension 1 x Iq, namely V^(/) = \V x,s{f)]s<^s(+)^ whose element 
yx,s{f) can be obtained from (see also (45a)) 

VxAt) = VtM W^l°' n {J27ThrrN.^L+^ S^'^ <^(t + (L - z)T)) ■ 

(75a) 



Nc-l 



Yl exp (j2TT hnN,+i y^'V(^ - 



i=0 



through an ordinary Fourier transform, that is Vx,s{t) — > Vx,s{f)- 

So, in the general case of non-integer modulation factors (p > 1) from (68) and (74) we obtain 

(r\ r / °° \ / n trr. ^ \ -1 1 (76) 

K\f)=T,Ko{f) + 2T,^ - - - - - 

where 



ij,-n){e='-f^^ij,-n{±) + n) k,U) 



Koif) = E P K = ^] Vxif) diag(p) Vx*{f) 

X 

K,{f ) = Ci(+) $*(/) = Y,V[Xrn = x] Ex{+) diag(p) K*(/) (76a) 

X 

K,{f) = $(/) C2(+) = E P K = ^] VM) 
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In the very special case of integer modulation factors {p = 1), we recall (67) to write the 
spectral lines as 

R^fif) = E \K,{kF,)pmf - kF,) (77) 

k 

while the continuous part of the spectrum in (76) needs a correction factor 

AR'^::^{f) = -\K,{f)W (78) 

according to (74). 



VI. Examples 

In this section we give some examples of PSDs evaluated through (76). We begin by showing 
some results for full response CPM signals (i.e. with L = 1) in Fig. 8. All plots refer to a 
continuous phase frequency shift keying (CPFSK) employing a phase response 

ro ,t<0 



t 



,0<t<LT (79) 
a ,t>LT 



and for different choices of the multi-/i sequence factors. Examples were taken from [9] and 
[10]. 

A further example of full response signaling taken from [9] is shown in Fig. 9, where CPFSK 
is compared to raised cosine (RC) signaling 

{0 ,t<0 
2IT - ^ sin (27r ^) ,0<t<Lr (80) 
I ,t>LT 

A final example is given for the partial response multi-/i CPM signaling formats of [12] in 
Fig. 10. The figure shows two RC formats, and one binary Gaussian minimum shift keying 
(OMSK) format with L = 4, having a phase response 

1 

^71 



^{t) = ^(k{^- D) -^(k{^- D) , I)(a) = a $(a) + ^e"^* 



where K = vr/ (2Vln2), and $(a) = -4= e 2* dt is the Gaussian normalized cumulative 



/2n 

distribution function. 
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VII. Conclusions 
Acknowledgements 
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Fig. 8 — Normalised PSD plots of full response CPFSK, for different multi-/i sequences, and for equally likely 
input symbols. 
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Fig. 10 — Normalised PSD plots for partial response single-/i and multi-/i CPM signals (equally likely input sym- 
bols). 
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Appendix 



A. Sketch of the proof of Proposition4 

We proceed step by step. For L = 1 we have 



_ Jjr„a 



which is a square matrix obtained by cyclically shifting the main diagonal to the left by r„a 
positions. The result can be obtained by evaluating the r„a power of the single step ciclical shift 
matrix = Dip- This corresponds to (34). 
For L = 2, we have 



where the first term can be further explicited as 



io,jo6A/'2p, iie.4]if 



jojo6A/'2p, jie.4jv/ 



\_^io,{ia+rn-il3)2p 



«o ,ioeA^2p 



® 



to obtain 

The prospect is similar for L = 3, where 



^«0,(j0+»'n-2il)2p 



«o,ioeA/'2p,iie^A 



«1 J2 



il,j2&AA 



-'12,0 



j26^A 



with the central matrix being an identity. The general result is thus 



Im 



L-2 



where the kronecker product of L — 2 occurrences of Im is simply I^l . As a consequence, (32) 
is valid. 
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B. Sketch of the proof of (36) 

Being 7r„ expressed as a Kronecker product in (33), we look for an eigenvector with the 
Kronecker structure 



2pxl Mxl Mxl 



providing the set of equations 



Wp (no ® ui) = no 



ImU2 = Ui 



q = Ul^i 

By solving the system we immediately obtain Ui 
equation 



n^_i = q, plus the remaining 



'2p Uo 



which is solved by no = ^ l2p, where l2p is a column vector of length 2p with all entries set to 
1. Note that the solution is valid independently on rn-L+i, and provides (36). 

C. Validity of (62) 

By exploiting the first of (45a) and (48), the matrix is 



Yx{i,j] 



,J65 



tV-^o r r. r . ... A 



a;jVc-i:*iVc+i:,-2 



io,.--,«JVc + I--2G-4i\/,ioGA/'2p,il,.--jL-ie-4ji 



where i = [io, . . . , 2Ar^+L_2]^ and j = [jo, • • • , Jl-i]^- It is now straightforward to see that the 
equivalent Kronecker form is given by (62). 

D. Sketch of the proof of PropositionQ 

For p > 1, we exploit the Kronecker products in (55) and (63) in (65), and observe that 
Vp Ip = 0. We then immediately obtain my{±) = thus assuring that no spectral lines are 
found, as stated in (66a). 
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Instead, for j» = 1 we have Vi li = 1 ■ 1 = 1, in which case it is easy to show that 
my(+) = q (g) ■ ■ ■ (g) q , my{-) = W2pmy{+) . 

Nc+L-l 

Then, from the equivalence 

ry^''\riT,) = myi+) my*(+) = my(-) my*(-) 
= [qq*] ® ■■■ ® [q q*] 



Nc+L-l 

we straightforwardly obtain (66b), by recalling (58). 
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